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ABSTRACT 

The jets of powerful blazars propagate within regions relatively dense of radiation produced 
externally to the jet. This radiation is a key ingredient to understand the origin of the high 
energy emission of blazars, from the X-ray to the 7-ray energy band. The main components 
contributing to the external radiation field are the accretion disk emission, including its X-ray 
corona, the broad line region, the infrared emitting torus and the cosmic background radia- 
tion. Their importance changes as a function of the distance from the black hole and of the 
value of the bulk Lorentz factor of the jet. These external radiation fields control the amount 
of the inverse Compton radiation with respect to the synchrotron flux. Therefore the predicted 
spectral energy distribution (SED) will depend on where the jet dissipates part of its energy to 
produce the observed radiation. We investigate in detail how the SED changes as a function 
of the location of the jet dissipation region, by assuming rather "standard" (i.e. "canonical") 
prescriptions for the accretion disk and its X-ray corona, the profile of the jet magnetic field 
and the external radiation. We confirm that most of the dissipation, if producing the 7-ray 
flux we see, must occur at hundreds of Schwarzschild radii from the black hole, to avoid the 
7-7 — ► process, and the consequent re-emission by the produced pairs. The magnetic en- 
ergy density of a "canonical" jet almost never dominates the radiative cooling of the emitting 
electrons, and consequently the inverse Compton flux almost always dominates the bolomet- 
ric output. This is more so for large black hole masses. Dissipation taking place beyond the 
broad line region is particularly interesting, since it accounts in a simple way for the largest 
inverse Compton to synchrotron flux ratios accompanied by an extremely hard X-ray spec- 
trum. Furthermore it makes the high power blazars at high redshift useful tools to study the 
optical to UV cosmic backgrounds. 
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1 INTRODUCTION 

Relativistic jets in blazars transport energy in the form of bulk 
motion of protons, leptons and magnetic field. When part of this 
power is dissipated, the particles emit the beamed radiation we ob- 
serve, consisting of two broad humps. The origin of the low fre- 
quency hump is well established, believed to be synchrotron radi- 
ation from relativistic (in the comoving frame) leptons. The nature 
of the high energy hump is a controversial issue. In general, we 
can group the existing models into three families: i) the high en- 
ergy radiation is generated by the same leptons producing the syn- 
chrotron, through the inverse Compton process (Maraschi, Ghis- 
ellini & Celotti 1992; Dermer & Schlickeiser 1993; Sikora, Begel- 
man & Rees 1994; Ghisellini & Madau 1996; Bloom & Marscher 
1996; Celotti & Ghisellini 2008); ii) There are two populations of 
leptons, one accelerated directly by the acceleration mechanism 
(i.e. shocks), and a second one resulting from cascades initiated 
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by ultra-relativistic protons (Mannheim 1993; Miicke et al. 2003; 
Bottcher 2007); iii) Ultrarelativistic protons emit by the proton- 
synchrotron process at high energies (Aharonian 2000; Miicke & 
Protheroe 2001). 

Another issue of debate is the role of electron-positron pairs. 
Sikora & Madejski (2000) discussed this problem concluding that, 
though they can exist, their amount is limited to a few pairs per 
proton, and a similar conclusion was reached by Celotti & Ghis- 
ellini (2008) analysing a large number spectral energy distributions 
(SEDs) for deriving the jet powers of blazars. This contrasts the 
idea, originally put forward by Blandford & Levinson (1995) that 
the 7-rays spectrum is the superposition of the spectra originating 
at different distances from the black hole, each one cutted-off at 
a different 7-ray energy as a result of photon-photon absorption. 
In this scheme a large amount of pairs is produced, since the high 
energy hump carries most of the emitted power: if this is partly ab- 
sorbed, we transform most of the total power into (energetic) pairs. 

Ghisellini & Madau (1996) argued that the Blandford & 
Levinson idea has one important observational consequence in 
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powerful blazars with a standard accretion disk: since the pairs are 
born relativistic, they contribute to the emission mainly by inverse 
Compton scattering the dense UV radiation field coming from the 
accretion disk. In fact, if we want the emission region be compact, 
it is reasonable to locate it close to the accretion disk and its X-ray 
corona. This implies that the coronal X-rays are targets for the 7- 
7 — > e process, and that the accretion disk UV photons become 
seeds for the scattering made by the newly born pairs. The result- 
ing radiation is mainly in the X-ray band, that should have a power 
comparable to the power absorbed in the 7-ray band, contrary to 
what observed. This led Ghisellini & Madau (1996) to conclude 
that the dissipation region in blazar jets cannot be very close to 
the accretion disk. On the other hand the observed fast variability 
argues for not too large distances. Taken together, these two lim- 
its strongly suggest that there is a preferred distance at which the 
'y-ray radiation is produced, at some hundreds of Schwarzschild 
radii. Since this argument makes use of the reprocessed radiation 
produced by the pairs, it cannot be applied when most of the radia- 
tion is produced below the pair production energy threshold and/or 
if the accretion disk is radiatively inefficient: in both cases very 
few pairs are created. Indeed, Katarzynski & Ghisellini (2007) pro- 
posed a jet model in which a dissipation close to the accretion disk 
resulted in a low level 7-ray emission, with no contradiction with 
the existing X-ray data. 

These arguments, although correct, are qualitative, and in our 
opinion a detailed computation of the reprocessing due to pairs is 
not yet present in the literature. Therefore one of the aims of the 
present paper is to derive some limits on the location of the dis- 
sipation region in the jet of powerful blazars. We will do it in the 
framework of the "leptonic" class of models [family i) mentioned 
above], and we limit our analysis to blazars having a "standard" 
(Shakura & Sunyaev 1973) accretion disk. This implies to consider 
Flat Spectrum Radio Quasars (FSRQs, and their likely parent popu- 
lation, FR II radio-galaxies) and not BL Lac objects, that are likely 
to have radiatively inefficient accretion flows (see Celotti & Ghis- 
ellini 2008 and Ghisellini & Tavecchio 2008 for more discussion 
about this point). 

While doing this, we will study the relative importance of dif- 
ferent sources of seed photons as a function of the distance Rdi BB 
of the dissipation region from the black hole, and their spectrum as 
seen in the comoving frame. This is the second aim of the paper. 
Besides the radiation coming directly from the accretion disk and 
its X-ray corona, we will consider the radiation produced in the 
Broad Line Region (BLR) and in a relatively more distant dusty 
torus, intercepting a fraction of the disk radiation and re-emitting 
it in the infrared. Finally, we also include the contribution of the 
Cosmic Microwave Background (CMB), important for very large 
(beyond 1 kpc) i?di ss - 

The term "canonical high power blazars" refers to the rather 
standard choice for both the environment of these sources and their 
jets: 



power, which is then conserved. After the acceleration phase, pos- 
sibly magnetic in origin, we assume that also the Poynting flux is 
conserved (see e.g. Celotti & Ghisellini 2008). 

Our study is not completely new, since several of its "ingredients" 
have already been discussed in the literature. The paper more ger- 
man to our is Dermer et al. (2009), but we include some new in- 
gredients. The novel features of our investigation concern mainly: 
i) the inclusion of the X-ray corona as an important producer of 
target photons for the 7-7 — > e process; ii) the calculation of 
the emitting particle distribution, including pair creation; iii) the 
effects of jet acceleration at small distances from the black hole; 
iv) the strict link between the properties of the accretion disk and 
the amount of the external radiation and v) the overall scenario al- 
lowing to describe in a more general way (than done before) the 
SED properties of high power jets at all scales. 

The paper is divided into four parts. In the first we study the 
different sources of external radiation, and the corresponding en- 
ergy densities, as seen in the comoving frame, as a function of dis- 
tance from the black hole and as a function of the bulk Lorentz fac- 
tor of the jet. In the second part of the paper we investigate the role 
of pair production processes when the dissipation region is close 
to the black hole, with the aim to find quantitative constraints. In 
the third part we construct the expected SED as a function of dis- 
tance, highlightening what are the relevant external seed photons. 
Finally, in the fourth part, we apply some of the above results and 
considerations when modelling the SED of some blazars, used as 
illustrative examples, and we check if our scenario can reproduce 
the phenomenological blazar sequence (Fossati et al. 1998). 

We use a cosmology with ho = Ha = 0.7 and 57m = 0.3. 
We also use the notation Q = 10 x Qx in cgs units, unless noted 
otherwise. 



2 SETUP OF THE MODEL 

Our model is characterised by the following setup. The accre- 
tion disk extends from R ln = 3i?s to R out = 500i?s (-Rs 
is the Schwarzschild radius) and is producing a total luminosity 
Ld = r]Mc 2 , where M is the accretion rate and 77 is the accretion 
efficiency. Locally, its emission is black-body, with a temperature 
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(i) 



Below and above the accretion disk there is a hot corona, emit- 
ting UV and X-rays with a luminosity Lx = fxLd- For sim- 
plicity, the corona is assumed to be homogeneous between 3 and 
30 Schwarzschild radii. The spectrum is assumed to be a cut-off 
power law: Lx(y) oc v~ ax exp{—v/h> c ). 

The broad line region (BLR) is assumed to be a shell located 
at a distance 



(i) The accretion disk is a standard, "Shakura & Syunyaev" 
(1973) disk. 

(ii) Above this accretion disk, there is an X-ray corona, emitting 
a luminosity less than, but comparable to, the luminosity emitted by 
the accretion disk. 

(iii) the BLR and the IR torus are located at distances that scale 
as the square root of the disk luminosity (Bentz et al. 2006; Kaspi 
et al. 2007; Bentz et al. 2008; see the discussion about this point 
made in Ghisellini & Tavecchio 2008). 

(iv) The jet is assumed to dissipate only a fraction of its total 



i?BLR 



1 n!7 r V 2 

10 L d ' 45 cm 



(2) 



It reprocesses a fraction /blr of Ld in lines, especially the hydro- 
gen Lyman-a line, and continuum. Following Tavecchio & Ghis- 
ellini (2008), we assume that the spectral shape of the BLR ob- 
served in the comoving frame is a black-body peaking at a factor 
L times the (rest frame) frequency of the Lyman-a line. 

We also assume the presence of a torus (see Blazejowski et al. 
2000; Sikora et al. 2002), at a distance 



Rm = 2.5 x 10 18 L\% cm 



(3) 
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reprocessing a fraction /ir of the disk radiation in the infrared. 
Note that both Pblr and Rm. scale as the square root of Ld: this 
implies that, in the lab frame, the radiation energy densities of these 
two components are constant, as long as Pdiss is smaller than these 
two radii. We emphasise that our treatment of the torus emission 
is approximate: it is likely that the torus itself is a complex struc- 
ture, possibly clumpy (Nenkova et al. 2008) with a range of radii, 
extending also quite close to the black hole, where the temperature 
is just below dust sublimation (i.e. ~ 1500 K). Our approach fol- 
lows in part the results of Cleary et al. (2007), finding weak signs 
of hot dust emission in the studied spectra, and partly is dictated by 
simplicity. 

The emitting region is moving with a velocity /3c correspond- 
ing to a bulk Lorentz factor V. We call Pdiss the distance of the dis- 
sipation region from the black hole. We consider either a constant 
r or include an acceleration phase of the kind (see e.g. Komissarov 
et al. 2007; Vlahakis & Konigl 2004): 



r = min 



/ _R_\ 1/2 



(4) 



When the acceleration phase is taken into account, the jet is as- 
sumed to be parabolic in shape, becoming conical when V reaches 
its maximum value (see e.g. Vlahakis & Konigl 2004). Calling r 
the cross sectional radius of the jet, and R the distance from the 
black hole, we have 



<>R 



r = 
r = ip R 



.1/2 



r > r 

1 • max 



(5) 



where tjj is the semi-aperture angle of the jet in its conical part, and 
the constant <f) is fixed by assuming that at the start of the jet we 
have r = Ro = 3Ps- The parabolic and conical parts of the jet 
connect at R — SRsF 2 ^^. 

The total power carried by the jet, Pj, was assumed in Ghis- 
ellini & Tavecchio (2008) to be related to the mass accretion rate, 
i.e. Pj = rijMc 2 . Since Celotti & Ghisellini (2008) found that Pj 
is greater than the disk accretion luminosity, i.e. Pj>-£/d, we as- 
sumed (in Ghisellini & Tavecchio 2008) that r/j was greater than 
the corresponding efficiency of transforming the accretion rate M 
in disk luminosity. When fitting the data of specific sources, we do 
not specify a priori the total jet luminosity, which is instead a result 
of the modelling (once we assume how many protons there are for 
each emitting lepton). Therefore Pj is not an input parameter in our 
scheme, but it is a quantity derived a posteriori. What we will spec- 
ify, instead, is the power injected in the dissipation region, that is 
related to the power carried by the jet in relativistic electrons (and 
by cold protons, once we specify how many protons there are per 
emitting electron). 

We assume a value of the magnetic field in the dissipation 
region, that corresponds to a Poynting flux Pb. If the jet is mag- 
netically accelerated, then the initial (i.e. close to Po) Pb should 
be of the same order of Pj, becoming Pb = £BPj =const when 
r reaches its maximum value. To describe the profile of Pb we 
assume the following prescription: 



Pb = 7rr 2 r 2 d7 B = Pj 



1 - 



r/3 



(1- 6B) 



max /-'max 



(6) 



In this way Pb = Pj initially, becoming a constant fraction cb of 
Pj when the jet is conical. Here Ub = B 2 /(8tt) is the magnetic 
energy density. 

The energy distribution of the particles responsible for the 
emission is derived through the continuity equation, assuming a 



continuous injection of particles throughout the source lasting for 
a finite time. This time is the light crossing time t C ross = rdiss/c, 
where rdi ss is the size of the emitting blob, located at the distance 
Pdiss from black hole. We always calculate the particle distribution 
at this time. The reason for this approach is suggested by the fast 
variability shown by blazars, indicating that the release of energy 
is short and intermittent. Besides, we believe that this approach is 
the simplest that can nevertheless describe in some detail the par- 
ticle distribution. In fact, it allows to neglect: i) adiabatic losses 
(important after rdiss/c, which is also the time needed to double 
the radius); ii) particle escape (again important for times longer 
than rdiss/c) and iii) the changed conditions in the emitting region 
(since the source is travelling and expanding, the magnetic field 
changes). 

High energy particles can radiatively cool in a time shorter 
than tcross- Let us call 7 coo i the energy of those particles halving 
their energy in a time t CIoaa . Above 7 CO oi, and at t — t CIoaa , the 
particle energy distribution N (7, tcross) can be found by solving 



+ Q(7) + P(i) = 



(7) 



where 7 is the cooling rate of a particle of energy 7771c 2 , Q(j) is the 
source term (i.e. the injection of primary particles) assumed con- 
stant in time, and P(7) is the term corresponding to the electron- 
positron pairs that are produced in photon-photon collisions. The 
formal solution of Eq.[7]is 



7V( 7 ) = 



r»-[Q(7) + P( 7 )]d7 



7 



7 > 7cooi 



When electrons with 1 < 7 < 7 coo i do not cool in t cr , 
proximate the low energy part of N(<y) with 

N(j) ~ tcross [Q(j) + P(7)], 7 < 7cool 



(8) 



we ap- 



(9) 



Note that, within our assumptions, the particle distribution of Eq. 
[8] and Eq. [9] correspond to the maximum N(j). The injection of 
primary particles Q(-y) is a smoothly joining broken power law: 



0(7) = Qo 



(7/7b) 



1 + (7/7b)- sl+s 



(10) 



where 7b is a break energy. The pair injection term P{"y) corre- 
sponds to the 7-7 —* process only, and it is calculated with the 
prescriptions given by Svensson (1987) and Ghisellini (1989). 

The total power injected in the form of relativistic electrons, 
calculated in the comoving frame, is 



Pi = m e c 2 V 



7Q(7)^7 



(11) 



where V = (47r/3)rdi ss is the emitting volume. Note that this is 
not equivalent to the power that the jet transports in the form of 
relativistic particles (as measured in the comoving frame), since P/ 
includes also the energy that will be emitted, and possibly trans- 
formed into pairs. 



3 ENERGY DENSITIES 
3.1 Direct disk radiation 

Each annulus of the accretion disk is characterised by a different 
temperature and it is seen under a different angle £ (with respect 
to the jet axis), thus its radiation is boosted in a different way. A 
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where Rx is the extension of the X-ray corona. 



3.3 BLR radiation 




Figure 1. Cartoon illustrating the accretion disk, its X-ray corona, the broad 
line region and a schematic representation of the IR torus. At the distance 
-Rdiss the jet is assumed to dissipate. At this distance, here assumed to be 
outside the BLR, we label the relevant angles for calculating the contribu- 
tion of the BLR radiation to the corresponding energy density. 



Within Rblr the corresponding energy density seen in the comov- 
ing frame can be approximated as (Ghisellini & Madau 1996): 



17F 2 /BLR,£d 



12 



Rdiss < Rblr 



(19) 



At distances much larger than -Rblr, and calling /i = cos a, we 
have (see the cartoon in Fig.[TJ 



4ttR 2 lrC 3j3 
(l-« 3 )] 



[2(1 - - (1 - /3 M2 ) 3 



/'i 

/I2 



— [1 + -Rblr/ Rdis 

= [i — -rIlr/RcL 



Rdiss 2> Rblr 

j-l/2 
,1/2 



(20) 



For i?diss^i?BLR the exact value of f/ B LR depends on the width 
of the BLR, which is poorly known. For this reason, in the range 
Rblr < Rdiss < 3Rblr we use a simple (power-law) interpola- 
tion. 

The BLR is assumed to "reflect" (Compton scatter) a frac- 
tion /blr.x (of the order of 1 per cent) of the corona emission. 
The existence of this diffuse X-ray radiation is a natural outcome 
of photo-ionisation models for the BLR (see e.g. Tavecchio & 
Ghisellini 2008, Tavecchio & Mazin 2009). The assumed value, 
/blr,x ~ 0.01, is the average value found for typical parameters 
of the clouds. 



stationary observer with respect to the black hole (lab frame) will 
see a flux, integrated over all annuli, given by 

2hv 3 /c* 



F d {v) = 2tt f I{y)dn = 2n f 



exp[hv/(kT)] 



-dfi 



(12) 



where fi = cos f , and ^d is given by 

Md = [i + RoWRLsr 1/2 d3) 

In the comoving frame of the blob, frequency are transformed as: 

v =bu, b = T(l - Pfi) (14) 
and solid angles transform as: 

dO dn 

(15) 



<m> = f = 2 .f 



where primed quantities are in the comoving frame. The intensities 
as seen in the comoving frame transform as: 

I' d (u') = b s I d {u) = b 3 I d (u'/b) (16) 

The specific radiation energy density seen in the comoving frame 
is 

= \f Uv'W = ^ f ^d/x (17) 



3.2 Radiation from the X-ray corona 

According to our assumptions, the total radiation energy density 
U' x of this component is (see e.g. Ghisellini & Madau 1996): 



U' x = 



fix 



f x L d r 2 

[1 + Rx/Rdis 



i-^x-p{i-A) + ^{i-A) 



-1/2 



(18) 



3.4 Radiation from the IR torus 

This component scales as [/blr, but substituting Rblr with Rir. 
We have 



17™ 



fm.L d r 
47TR?,, c 



Rdiss < Rir 



(21) 



For Rdiss > Rir we have the same behaviour as in Eq. [20] but 
with Rir replacing Rblr. In a vF v plot the (lab frame) peak fre- 
quency of this component is assumed to be at i^ir = 3 x 10 13 
Hz (see Cleary et al. 2007), independent of the disk luminos- 

1/2 

ity, since Rir scales as L d . The corresponding temperature is 
Tjr = ft^iR/(3.93fc) (we must use the factor 3.93, instead of the 
usual 2.82, because we are using the peak frequency in vF v ). In the 
comoving frame this corresponds to 



370 b K 



(22) 



3.5 Radiation from the host galaxy bulge 

The bulge of the galaxy hosting the blazar can be a non-negligible 
emitter of ambient optical radiation (see e.g. (Stawarz, Sikora & 
Ostrowski 2003). Within the bulge radius R s tar emitting a lumi- 
nosity Lstar we have 

t^2 -^star 



UL 



4ttR 2 



(23) 



As an order of magnitude estimate, we have (7 s ' tar ~ 10~ 10 r 2 erg 
cm' 3 using Rstar = 3 x 10 44 erg s _1 produced within a bulge 
radius R s tar ~1 kpc. When Rdiss > Rstar, t/star decreases in an 
analogous way as (7/ R and )7blr (once we substitute Rstar to Rir 
or Rblr). 
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3.6 Radiation from the cosmic background 

The energy density of the Cosmic Background Radiation (CMB), 
as seen in the comoving frame, is 



^CMB = aT 4 r 2 (l + zf 



where a = 7.65 x 10 15 erg cm 3 • 
temperature of the CMB now (i.e. at z — 0). 



(24) 

4 and To = 2.7 K is the 



3.7 Magnetic field 

Following our prescriptions concerning the power carried in the 
form of Poynting flux, we have 



Vk = 



T(3 



,i3n 



maxpmax 



(25) 

The magnetic field B scales as 1 /-Rdiss both in the acceleration and 
in the coasting (i.e. V =const) phases, but with a different normal- 
isation. 



3.8 Internal radiation 

This component corresponds mainly to the radiation produced by 
the blob that can be efficiently scattered through the inverse Comp- 
ton process. This cannot be calculated without specifying the rela- 
tivistic particle distribution. If the emitting volume V is a sphere, 
the synchrotron radiation energy density formally is: 



U' 

w syn 



V 



4 r 1rr , JjV(7)7 2 d 7 

-[nr diss a T ]U B f 



-[nr diBB a T (-y 2 )]UB 



(26) 



where err is the Thomson cross section, n is the number density 
of the emitting particles and y = (XTWdiss^ 2 ) is the relativistic 
Comptonization parameter. This parametrisation follows the one in 
Ghisellini & Tavecchio (2008). From the fits to a sample of blazars 
performed in Celotti & Ghisellini (2008) we have that y is in the 
range 0.1-10. Since the found distribution of y values is rather nar- 
row, one has a first order estimate of what is the internal radiation 
through the value and profile of the magnetic energy density Ub ■ 



3.9 Comparing the different components 

We have calculated the contributions to the energy density given 
by the different components as a function of Rdi ss < plotting them 
in Fig. [2] For these cases, we have chosen a black hole mass 
M = 1O 9 M0, an accretion disk extending from 3R B to 500 Rs, 
an accretion efficiency r] — 0.08, /blr = fx = 0.1, fm = 0.5. 
We plot the results for two values of L d ' 10 45 and 10 47 erg s _1 . 
In the top panel we show the case of an accelerating jet, whose 
bulk Lorentz factor is F — min[15, (Rdi ae /3Rs) 1 ' 2 ]. In the bot- 
tom panel we assume a constant T = 15 along the entire jet. As 
can be seen, the dominant energy density is different at different 
-Rdiss- As a rule, all the external radiation energy densities drop 
when i?diss is greater than the corresponding typical size of the 
structure producing the seed photons. 

Particularly interesting is the comparison between U d and 
[/blr- The distance above which £/ BLR > U d depends upon the 



R di SS / R S 
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6 
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Rdiss [cm] 



10 20 



Figure 2. Top panel: Comparison of different energy densities as mea- 
sured in the comoving frame. The moving blob is assumed to have a bulk 
Lorentz factor V = min [15, (i?di ss /3R s ) 1/2 ]- The black hole has a mass 
M = 1O 9 M0. The different contributions are labelled. The disk emits 
as a blackbody, and extends from 3 to 500 Schwarzschild radii. The two 
sets of lines correspond to two disk luminosities: 10 45 and 10 47 erg s — 1 . 
The radius of the broad line region is assumed to scale with the disk lumi- 
nosity as -Rblr, = 10 17 i d ^45 C1TL T fie ra ^ ms °f tne torus scales as 
Rjji = 2.5 X 10 18 Lj //2 g cm. The X-ray corona is assumed to be homo- 
geneous, to extend up to 30 Schwarzschild radii and to emit 10 per cent of 
the disk luminosity. The contribution of the BLR between 1 and 3 -Rblr. 
depends on the unknown width of the BLR itself (dotted line). The mag- 
netic energy density (long dashed lines) is calculated assuming P; = L d 
and ee = 0.1. Bottom panel: as above, but assuming a constant T = 15 all 
along the jet. 
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10 



R di SS / R S 

100 1000 



10 4 



10 5 



m 




R diss [cm] 

Figure 3. The energy density as seen in the comoving frame of the blob, 
moving with a bulk Lorentz factor T = min[15, (iJdiss/3-Rs) 1 ^ 2 ]- The 
black hole has a mass M = 10 9 A/q, and the disk emit -Ldisk = 10 46 erg 
s _1 . The different contributions are labelled. This figure shows the effect 
of changing the outer radius of the accretion disk, from 150 to 5000 R<j, 
as indicated by the labels and arrows. As the outer radius increases, U' d 
increases at distances -Rdiss ~ Rout, but there the disk contribution is 
overtaken by l/ BLR and t/j' R . 



disk luminosity, since the BLR adjusts its radius so to give a con- 
stant (in the lab frame) energy density. So for L d = 10 47 erg s _1 
the energy density from the BLR dominates only above 10 17 cm, 
equivalent to ~ 300-Rs, while for Ld = 10 45 erg s _1 it starts 
to dominate about three times closer, when V has not reached yet 
its maximum value. Note that, for the shown cases, dominates 
over the external radiation energy density only at the start of the jet 
and up to 100-300 Rs, where E/ B lr takes over. In this particular 
examples, we assumed Pb = O.lLd- 

Fig.[3]illustrates the effects of changing the outer radius of the 
accretion disk, from -R ou t = 150 to 5000 Schwarzschild radii. As 
expected, U' d remains the same at small and very high distances, 
and increases with Rout in between, at distances comparable with 
Rout - The increase of U' d is modest, and occurs when the external 
radiation energy density is dominated by the BLR and torus com- 
ponents. 

In Fig. [4] we show U' d as a function of V for different loca- 
tions of the dissipation region. Of course the largest U' d are obtained 
when the source is very close to the accretion disk. For each i?diss 
there is a value of F that minimises U' d . To understand this (some- 
what anti-intuitive) behaviour we must recall that the boosting fac- 
tor b = T(l — /3/i) has a minimum for /3 = fi. This corresponds, 
in a comoving frame, to // = 0, i.e. an angle of 90°. Therefore, 
when the external radiation is all coming from a ring, we expect 
that in the comoving frame, the radiation energy density initially 
decreases increasing 0, but when the photons are starting to come 
from the forward hemisphere (i.e. // < 0, or, equivalently, fi < /?), 
U' d increases increasing (3. This explains the minima shown in Fig. 




1 10 100 

r 



Figure 4. The radiation energy density U' d (as seen in the comoving frame) 
produced by the direct radiation of the accretion disk as a function of T, for 
different distances -Rdiss. f° r a total disk luminosity L d = fO 46 erg s — 1 
and a black hole mass M = 10 9 Mq. For comparison we show the radi- 
ation energy density observed within the BLR. It can be directly compared 
with U' d only for _R d i ss < Rblr- 



[4] once we weight the above effect with the intensity produced at 
each ring. 



3.10 Spectra of the external radiation 

Besides comparing the frequency integrated components forming 
the radiation energy density in the comoving frame, it is instruc- 
tive to compare also their different spectra. This is because the 
Doppler boost appropriate for each component is not the same, 
and it changes for different distances from the black hole. In Fig. 
[5] we show one illustrative example, calculated for a blob mov- 
ing with a constant V = 15 and located at different 7?diss- The 
disk is assumed to emit L d = 10 45 erg s _1 and the black hole 
mass is M = lO 9 A/0, corresponding to Rs = 3 x 10 14 cm. 
With this disk luminosity, the BLR is located at -Rblr = 10 17 
cm, implying that for -Rdiss < Rblk the radiation energy den- 
sity from the BLR is the same. For larger -Rdiss the total t/ BLR 
decreases (see Eq.|20t. The radiation energy density from the disk 
and its corona always decreases increasing -Rdiss, being larger than 
[/ BLH for i?diss < 100i?s ~ 3 x 10 le cm (see also Fig. HJ. 
The disk radiation dominates again for R dlss > x 10 18 cm (see 
Fig. |2). For small -Rdiss, the high energy X-ray photons com- 
ing from the corona will be responsible of a large optical depth 
for the 7-7 process, as discussed below. The IR energy density 
coming form the torus is constant for all Rdiss but for the largest 
value (since R m = 8.3 x 10 3 7?s = 2.5 x 10 18 cm). This com- 
ponent will be larger than U' d and (7 B lr f° r Pdiss greater than 
~ 500R S ~2x 10 17 cm. 
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Figure 5. The spectra of the radiation energy density (as seen in the co- 
moving frame) produced by the direct disk plus the X-ray corona (t7j, 
solid line), the BLR (^b LR , including the "reflected" X-rays, dashed line) 
and the IR torus (C/j' R , dotted line). We have assumed a constant T = 15, 
M = 10 9 A/q and = 10 45 erg s — 1 . The spectra are shown for 5 dif- 
ferent Rdiss- F° r Rdiss = 10 an d 100 Schwarzschild radii U^„(v') is 
the same, since for these cases -Rdiss < Rblr- The contribution of the 
IR emitting torus is the same up to -Rdiss = 10 3 i?s, and starts to decline 
afterwards, when -Rdiss > Rir- 

3.11 Input parameters 

Having included many components of external radiation, and hav- 
ing tried to link the jet emission to the accretion disk luminosity and 
the black hole mass, we necessarily have a large numbers of param- 
eters for a specific model. It can be useful to summarise them here, 
dividing them into two separated lists. In the latter we have those 
parameters that we fix for all models, based on physical considera- 
tions or a priori knowledge, or that do not influence the model SED 
(apart from pathological cases); in the former we list the important 
input parameters. 

• M: the mass of the black hole; 

• Ld : the luminosity of the accretion disk. 

• R<Ubb'- the distance from the black hole where the jet dissi- 
pates. 

• P{ : the power, calculated in the jet rest frame, injected into the 
source in the form of relativistic electrons. 

• B: the magnetic field on the dissipation region; 

• r max : the value of bulk Lorentz factor after the acceleration 
phase, in which the bulk Lorentz factor increases as [R/^Rs)] 1 ^ 2 ■ 

• 9 V : the viewing angle. 

• 7b: value of the random Lorentz factor at which the injected 
particle distribution Q(y) changes slope. 

• 7m ax : value of maximum random Lorentz factor of the in- 
jected electrons. 

• S2' slope of Q("/) above 7b. 

These amount to 10 parameters. If we use V ~ 1/r, as dis- 
cussed in Sj9] then the number of relevant parameters decreases to 
9. Note that we do not consider Pj as an input parameter, since it 



is derived once calculating the total electron energy density in the 
dissipation region, and assuming a given number of (cold) proton 
per emitting electron. 

The following parameters are either unimportant or are well 
constrained: 

• Rx : the extension of the X-ray corona. We fix it to 30-Rs- 

• -Rout: the outer radius of the accretion disk. For the shown 
illustrative examples, we fixed it at 500 Schwarzschild radii, and 
we have shown that changing it has a modest influence on U' d . 

• s\\ slope of Q(y) below 7b. It is bound to be very flat (i.e. 

-1< 81 < 1). 

• ax : the spectral shape of the coronal X-ray flux. We fix it to 

ax = 1. 

• hv c : the high energy cut-off of the X-ray coronal flux. We fix 
it at 150keV. 

• Lx : total X-ray luminosity of the corona. We fix it to Lx = 
0.3Ld. 

• /blr: fraction of Ld intercepted by the BLR and re-emitted 
in broad lines. We fix it to /blr = 0.1 

• /blr ,x: fraction of Lx of the corona scattered by the BLR. 
We fix it to /blr,x = 0.01 

• /ir : fraction of Ld intercepted by the torus and re-emitted in 
IR. We fix it to /ir = 0.5. 

• tp: semi-aperture angle of the jet. We fix it to ip = 0.1. 

• C/star: this is the radiation energy density within the bulge 
of the host galaxy. We approximated it with the constant value 
[/star = 10" 10 erg cm" 3 . 



4 THE ROLE OF PAIRS 

At large distances from the black hole, the density of the photon 
targets for the 7-7 — > process is small, with few pairs being 
produced inside the emitting region. At small distances, instead, the 
presence of the X-ray corona makes pairs easy to create. With our 
assumptions, we can see the effects of the reprocessing of the spec- 
trum due to pairs. Whenever pairs are produced inside the emitting 
region at small -Rdiss, the radiative cooling rates are large, implying 
that all particles cool in one light crossing time. This in turn implies 
a particle distribution of the form given by Eq.[8] with 7 coo i ~ 1. 

One necessary ingredient for making pairs important is the rel- 
ative amount of power emitted above the energy threshold m e c 2 (in 
the comoving frame). If this power is a small fraction of the total, 
the reprocessing will be modest, even if all the 7-rays get absorbed. 
There is then the possibility to have a large jet power dissipated 
close to the accretion disk, but with electrons of lower energies, 
whose emission is mostly below threshold. Furthermore, even if the 
electron energies are large, the bulk Lorentz factor could be small, 
especially very close to the accretion disk, where presumably the 
jet just starts to accelerate. In this case the magnetic energy density 
becomes more dominant with respect to the disk radiation, imply- 
ing a relatively modest high energy emission, and furthermore the 
whole spectrum is much less boosted. 

In the following we will consider the case of two jets with 
the same parameters, except that the first is moving with a large T 
already at small J?diss, the other, instead, is accelerating. 

Consider Fig. [6] it shows the effects of including the 7- 
7 — > e on the final spectrum. In the middle panel of Fig. [6] the 
same spectra are shown as seen in the comoving frame. In this 
case on the y-axis we plot the specific compactness, defined as 
t' (v') = ctL' {v') / {r d i SB m e (?). The bottom panel of this fig- 
ure shows the particle distribution (considering or not the produced 
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Figure 6. Top panel: example of how reprocessing due to electron-positron 
pairs can be important, if the dissipation region is located too close to the ac- 
cretion disk and its X-ray corona, and if this region produces 7-rays above 
the pair-production threshold. The two humps at low energies correspond 
to the flux produced by the accretion disk and the IR flux from the torus. 
See Tab.[T]for the parameters used to construct the shown SED. Mid panel: 
the SED as observed in the comoving frame. Bottom panel: the particle dis- 
tribution 7 3 t(7). For this example, the maximum particle Lorentz factor is 
7max = 10 4 and 7t> = 10 3 . A redshift 2 = 3 has been assumed. 



pairs) in the form 7 3 t(7) = 7 3 o-Tf"dissA r (7). The 7 s factor makes 
the 7 3 t(7) distribution to have a peak: electrons at this peak are 
the ones producing the peaks in the vL v synchrotron and inverse 
Compton spectra. The parameters used for this example are listed 
inTab.Cn 

One can see that pairs redistribute the power from high to 
low frequencies. The N(-/) distribution steepens at low energies 
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Figure 7. The compactness £' = 47rf/'o"Tfdiss/( rn eC 2 ) as measured in 
the comoving frame of a source of size r<jj aa as a function of i?diss ■ We use 
the same parameters used to construct the top panel of Fig. [2] 



(where pairs contribute the most), and this softens the emitted spec- 
trum especially in the X-ray range. The reprocessing of the X-ray 
spectrum is stronger than at lower frequencies because the low en- 
ergy electrons emit, by synchrotron, self absorbed radiation, whose 
shape and normalisation are largely independent of the slope of 
N(-y). Instead, in the X-ray range, the presence of the intense 
optical-UV radiation field coming from the accretion disk implies 
that the dominant process is the external Compton (EC) scatter- 
ing with these photons. In turn, this implies that the radiation we 
see in the soft and hard X-ray range is due to relatively low en- 
ergy electrons. In fact, the scattered frequency v ~ 7 2 T 2 ^o, where 
vq is the (observed) frequency of the seed photon. Therefore, at 
10 keV, we are observing the radiation produced by electrons with 



7 ~ 50 !/ c 



-1/2 T 



' a few. 



In the above example we have, on purpose, assumed that the 
dissipation region, although very close to the black hole, is already 
moving at large speeds, corresponding to a bulk Lorentz factor of 
10. We may ask, instead, what happens if the same kind of dissi- 
pation occurs while the jet is still accelerating, and has therefore a 
relatively small bulk velocity when it is close to the black hole. 

To this aim, Fig.[7]shows the compactness I', as measured in 
the comoving frame, as a function of -Rdias, for an accelerating jet. 
We define this comoving compactness as: 

47T(7 / cr T rdi ss 



(27) 



where rdi ss is the size of the blob. The term U' is the external ra- 
diation energy densities. The compactness is directly associated to 
the optical depth of the 7-7 — > e^ 1 process, and it describes the 
absorption probability of a 7-ray photon while it is inside the blob. 
In other words, this definition does not account for the absorption a 
photon can suffer while it has already escaped the blob. It also ne- 
glects the contribution of the internal radiation energy density due 
to the synchrotron and the synchrotron self Compton (SSC) emis- 
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sion. Fig.Qshows that there are three distances where the compact- 
ness is large: 

(i) At the base of the jet the main targets for the 7-7 process are 
the X-ray photons produced by the corona. They can interact with 
7-ray photons just above threshold, implying that the absorption is 
important. However, at these distances, the magnetic energy density 
is larger than the external radiation energy density (dominated by 
the disk radiation), and this implies a modest external Compton 
scattering, and a modest production of 7-rays. 

(ii) The compactness becomes large again for Rdiss ~ Rblr- 
Since the main contribution to U' is due to the BLR photons, (seen 
at UV frequencies in the comoving frame), the absorbed 7-rays 
will have GeV energies (see next subsection for more details). 

(iii) The third relevant range of distances corresponds to 
Rdiss ~ Rir- The main targets for the absorption process are the 
IR photons produces by the torus, absorbing TeV 7-rays. 

Fig.[8]shows three SED corresponding to -Rdiss = 10, 100 and 
10 3 Rs, with an electron injection function equal to the one used in 
Fig. [6] but assuming V = min[15, (R di ss/3-Rs) 1/2 ]. Tab. H lists 
the parameters. We can see the dramatic changes with respect to 
Fig. [6] and also among the three shown models. Pair production, 
although present (compare solid and dashed lines), is marginal and 
there is no noticeable reprocessing of the primary spectrum. This 
is due to the much reduced importance of the external radiation 
when r is small. As a consequence, the spectrum is dominated by 
the synchrotron emission, with a small fraction of the power being 
emitted above the pair production energy threshold. This in turn 
implies that early (i.e. close to the black hole) dissipation is not 
an efficient mechanism to produce electron-positron pairs. Since 
a small V means small Doppler boosting, the produced spectra at 
small distances are undetectable, overwhelmed by the emission of 
more distant and fastly moving components. 

When r is large, therefore for the -Rdiss = 10 3 Rs case, the 
relevant absorbing photons are those of the broad line region, with 
a very modest contribution from the coronal X-rays reflected by 
the BLR material. The overall effect is marginal. 

We can conclude that a "canonical" jet can dissipate part of 
its kinetic energy even at small distances. One way to avoid strong 
pair production is to have small electron energies, emitting a small 
power above the pair production threshold. Alternatively, if the 
jet is accelerating, and close to the black hole it has a modest 
T-factor, the produced radiation is hardly observable because the 
Doppler boost is limited. In relative terms, the synchrotron lumi- 
nosity should be dominant, with a small fraction of power being 
emitted at high energies. 



4.1 Pair production versus Klein-Nishina effects 

In powerful blazars, the high energy flux, at ~GeV energies, is 
dominated by the inverse Compton process between high energy 
electrons and external radiation. In Tavecchio & Ghisellini (2008) 
we pointed out that when the dominant contribution to the seed 
photon for scattering is given by the BLR, the hydrogen Lyman-a 
photons are the most prominent ones. The fact that there is a char- 
acteristic frequency of the seed photons allows to easily calculate 
when the Klein-Nishina effects (decrease of the scattering cross 
section with energy) are important. For completeness, we briefly 
repeat here the argument, with the aim to extend it to cases in which 
Rdiss is beyond the BLR, and as a consequence, the relevant seed 
photons becomes the IR ones. 

In the comoving frame, and as long as we are within the 




15 20 
Log v [Hz] 

Figure 8. Same as Fig. [6] but for an accelerating jet and for different -Rdiss- 
In this case, in the comoving frame of the jet, the magnetic energy density 
is dominant for -Rdiss =10 an d 100Rg. Furthermore, the much reduced 
Doppler boosting for small values of Rdi ss implies a very small received 
flux. There is some 7-7 absorption (compare solid and dashed lines, the 
latter corresponding to switching off the 7-7 process), but involving a very 
small amount of power. Consequently, pairs reprocessing is unimportant. 
Light grey lines (green in the electronic version) correspond to the non- 
thermal spectrum, black lines include the contribution from the accretion 
disk, its X-ray corona and the emission from the IR torus. 

BLR, The observed seed Lyman-a frequency is observed at v^ a — 
2Yv\ ja . To be in the Thomson scattering regime we require that 



2Yhv^ a < 



7 



(28) 



If the random Lorentz factor 7 of the electron satisfies this condi- 
tion, the energy of the scattered photon is 

4 , 2Y8 _ S 



vkn = -7 — — 

3 1 + 2 



= 15: 



GeV 



(29) 



r(i + «) 

Above this energy, the spectrum steepens due to the decreased effi- 
ciency of the scattering process. This may well mimic the effect of 
photon-photon absorption, making it difficult to discriminate be- 
tween the two effects. Consider also that the Lyman-a photons are 
the best targets for the photon-photon process at these energies and 
in these conditions (namely, relatively far from the accretion disk, 
but still within the BLR). 

We can repeat the very same argument assuming that the rel- 
evant seed photons are the ones produced by the IR torus. This 
occurs when the dissipation region is beyond the BLR but within 
Rir. In this case the relevant frequency is v>m = 3 x 10 13 Hz and 
we can rewrite Eq.|29|as 



2YS 



"KN = -7 ^IR — 



1.2 



r(i + *) 



TeV 



(30) 



This implies that, if we do see a spectrum which is unbroken even 
above 15/(1 + 2) GeV, we can conclude that the dissipation region 
is beyond the BLR. 

Also in this case, the "intrinsic" absorption due to the 7-7 pro- 
cess is expected to become important at ~ ^kn- There is therefore 
a range, between 15/(1 + z) and ~ 1000/(1 + z) GeV, in which 
the intrinsic spectrum should not suffer from the Klein-Nishina nor 
from the 7-7 process due to "internal" absorption. These sources 
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are thus the best candidates for studying the photon-photon absorp- 
tion due to the cosmic optical-UV background radiation. 



5 DISSIPATION AT LARGE DISTANCES 

Although there are strong indications that most of the observed ra- 
diation originates in a well localised region of the jet, there is also 
compelling evidence that the jet dissipates also in other regions, 
larger and more distant from its base. In fact, in the radio band, 
we do see bright knots at different distances contributing to form 
the characteristic flat total radio spectrum of the entire jet, and we 
do observed optical and X-ray knots also at hundreds of kpc from 
the jet apex. The power associated to these, often resolved, features 
is a small fraction of the total observed jet bolometric luminosity, 
but in specific bands the flux produced at larger distances can be 
comparable to the flux produced in the most active region. 

It is therefore instructive to calculate the predicted spectrum as 
a function of distance, extending our analysis well beyond the re- 
gion of influence of the IR flux produced by the torus, and therefore 
analysing the effects of the radiation energy density of the cosmic 
background radiation (CMB). This has already been pointed out 
as an important source of seed photons for enhancing the inverse 
Compton flux of the bright X-ray knots detected by the Chandra 
satellite at tens and hundreds of kpc from the core of the jet (Tavec- 
chio et al. 2000; Celotti Ghisellini & Chiaberge 2001). 

Fig. [9] shows the different contributions to the radiation en- 
ergy densities as seen in the comoving frame, compared to U B , 
and including Uc MB , important for large distances. To construct 
this figure we have assumed M = 1O 9 M0, Ld = 10 47 erg 
s"\ P' = 10 44 erg s~\ Pb = 10 46 erg s"\ z = 3 and 
r = min[15, (-Rdiss/3-Rs) 1 ^ 2 ]- This figures shows that U B domi- 
nates only in two regions of the jet: in the vicinity of the black hole 
and at i? d iss ~ 10 20 -10 21 cm. Within these two distances, [/ BLH 
and Uxr dominate, and beyond 10 21 cm Uc MB takes over. This has 
two important and immediate consequences: 

(i) The external Compton radiation will be more important than 
the synchrotron (and the SSC, if y ~ 1) radiation at all distances, 
except the (relatively narrow) distance intervals where U B domi- 
nates. 

(ii) At large distances the jet is conical, and the blob size scales 
linearly with -Rdiss (namely the jet is a cone with the same aperture 
angle ip), Then the light crossing time t CIOSS — ipR^iss/c, This hy- 
pothesis means that 7 coo i reaches a maximum where U B = Uq MB . 

To understand the second issue, consider the value of 7 coo i after a 

time tcross • 



7cool — 



3m c c 



(3D 



In the regions of interest (i.e. above 1 kpc), Uq MB is constant, while 
U B oc rjjj^. Therefore the maximum 7 coo i occur when Uc MB = 
U B , i.e. at 



ipRc, 



Pb 



1/2 



T§(l + z) 2 F 2 



(32) 



Inserting this in Eq.[3T]and setting U' = Uc MB + U B = 2!7cmb> 



we obtain: 



max 
7cool 



3m c c 



4<7 T T 2 (l + z) 2 
3.08 x 10 6 



7TC 



1/2 



(1 + Z) 2 Pj 



1/2 
B,46 



(33) 



Independent of V. Electrons with this energy emit an observed syn- 
chrotron frequency ^ CO oi given by: 



3.6 x 10 6 B 7c 2 ool 5 



1.09 x 10 1 



r<5 



(l + z) 3 L B , 



Hz 



(34) 



Analogously, the observed frequency emitted by these electrons 
scattering the peak of the CMB radiation is 



ic 

^cool 



-7coolI 



2.8 x 10 2 



T6 



(1 + z) 4 P B ,46 



Hz 



(35) 



We thus expect that the synchrotron spectrum produced in pow- 
erful jets, at large distances, cuts-off at a frequency given by Eq. 
1341 or somewhat larger if, at these scales, the active region has a 
size smaller than ipR (implying a smaller cooling time, and thus 
electron energies greater than the ones given by Eq.l33t. 

Fig. 1 101 shows the predicted SEDs corresponding to different 
-Rdiss (indicated by the vertical grey lines in Fig. [9] with the corre- 
sponding numbers); all the relevant input parameters for the differ- 
ent SEDs are reported in Tab. Q] The profile of the external radia- 
tion and the value of the magnetic field correspond to what shown 
in Fig. [9] We have also assumed that, increasing Rdiss, the injected 
power in relativistic electrons decreases by a factor 3 increasing 
Rdiss by a factor 10. At the same time, we assumed that the max- 
imum energy of the injected electrons increases by a factor 3 for 
a tenfold increase of -Rdiss- These choices are arbitrary, but reflect 
the observational evidence that the bolometric radiative output of 
jets decreases with distance. Furthermore the existence of optical 
jets, whose emission is due to the synchrotron process, ensures that 
at large jet scales there are very energetic electrons. 

The first three SEDs (number 1, 2 and 3, corresponding to 
Rdiss = 10 3 , 10 4 and 10 5 Rs, respectively) have a large Compton 
dominance, corresponding to the large ratio between the external 
radiation and the magnetic energy densities. Note that only the first 
uses the BLR photons as the main seeds for the Compton process 
(and thus it has the largest peak frequency), while the SED 2 and 
SED 3 use the IR photons as seeds. SED 1 is dominating in the 7- 
ray band (above ~ 1 MeV), while in the far IR to UV bands there 
is the contribution of the thermal radiation (accretion disk and IR 
torus). Remarkably, the soft X-rays are produced almost equally 
by these three jet dissipation sites (see also below). 

SED 4 (i.e. Rdiss = 10 6 Rs) is the only one where the mag- 
netic energy density dominates. Correspondingly, the synchrotron 
flux dominates the bolometric output. Comparing SED 3 and SED 
4, we note that they have total luminosities that differ by more than 
a factor 3. This is due to the incomplete electron cooling, and im- 
plies that not all the power injected in random energy of the elec- 
trons can be radiated in one crossing time. For the same reason, 
there is a rather large jump in 7 coo i between -Rdiss = 10 and 
10 6 i?s (see the bottom panel of Fig. [9), corresponding to the fast 
drop of the energy density of the external radiation between these 
two distances. 

SED 5 (-Rdiss = 10 7 i?s, or 1 kpc) has again a relatively large 
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Figure 9. Top panel: the contributions to the radiation energy density as 
seen in the comoving frame of the emitting blob, as labelled, and the mag- 
netic energy density as a function of the distance of the blob to the black 
hole. We assumed: M = 10 9 M Q , L d = 10 47 erg s" 1 , z = 3 and 
T = min[15, 0Rdiss/3Rs) 1/2 ]- The Poynting flux is P B = 10 46 erg 
s~ 4 once the jet has reached its maximum bulk Lorentz factor. The grey 
vertical lines indicates the distances assumed to construct the SEDs in Fig. 
1 101 the different numbers help to identify the corresponding SED. The bot- 
tom bottom panel shows 7 coo i (after one light crossing time since the start 
of the injection) as a function of i?diss- 

Compton dominance, due to the prevailing of the CMB and the 
starlight (from the galaxy bulge) energy densities over the magnetic 
one. This SED has almost the largest possible 7 coo i (see Fig. [9} 
and this is reflected in the calculated SED, having the largest peak 
frequency of the high energy component. 

SED 6 and 7 (i? d i ss = 10 8 and 10 9 -R S , or 10 and 100 kpc, 
respectively) have a very large Compton dominance, due to the fact 
that Uq MB is constant, while Ub decreases with distance. As dis- 
cussed previously, the constancy of Uc MB makes 7 CO oi to decrease 
with -Rdiss, since the cooling time at which 7 coo i is calculated in- 
creases. This also implies that, despite the decrease of P{, the SED 
5, 6 and 7 have the same bolometric luminosity, because a larger 
portion of the electron population can cool in one t cr oss- 

The sum of all SEDs is shown by the black thick line in Fig. 
[T7J1 Note that, in the radio band, the total flux F v oc v° (or slightly 
harder), as observed. At frequencies greater than the radio ones, the 
total flux is dominated by SED 1, but the flux originating at rela- 
tively large scales can be important in some frequency bands, as 
the soft X-ray one, with important consequences on the observed 
variability and the correlation of variability in different bands. Ac- 
cording to Fig.[l0](which is, we re-iterate, only one possible exam- 
ple, shown for illustration) the correlation between the 7-ray and 
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Figure 10. Sequence of SEDs calculated for different i?diss from 10 3 i?s 
to 10 9 Rg (one per decade). The injected electron luminosity is P. = 10 44 
erg s — 1 for -Rdiss = 10 3 Rs an d is reduced by a factor 3 each decade. The 
particle distribution has always the same 7j, = 100, while 7 ma x increases 
by a factor 3 each decade starting from 7 ma x = 10 4 for iJdiss = W 3 Rs- 
All other parameters are the same as in Fig. [9] What shown are the observed 
spectra neglecting the absorption of the high energy flux due to the IR- 
opt-UV cosmic background. The numbers correspond to the same numbers 
in Fig.|9]and correspond to the SED at different -Rdiss- The thicker black 
line is the sum of all the SED. The grey line at radio frequencies indicates 
F(u) oc u° . The received flux is calculated assuming that the source is at 
z = 3. 

the IR flux should be tighter than the correlation between the 7- 
ray and the soft X-ray flux, diluted by the the flux originating at 
larger scales (that can vary on longer timescales). Furthermore, if 
the soft X-ray flux is produced by external Compton, the energies 
of the electron emitting it are relatively modest (namely 7 ~1-10) 
while, in the 7-ray band, we see the emission of the electrons with 
the highest energies (which emit, by synchrotron, at IR-optical fre- 
quencies). The cooling timescales are therefore different, produc- 
ing different variability if i coo i at low energies is longer than the 
light crossing time. If not, then we expect the same variability pat- 
tern for emission produced in the same zone. One example of dif- 
ferent variability behaviour is displayed by the blazar 3C 454.3, as 
recently studied by Bonning et al. (2008), and has been interpreted 
by these authors as a consequence of the different energies of the 
electron contributing in the X-ray and 7-ray bands. The contribu- 
tion to the X-ray flux by other zones of the jet can be an alternative 
possibility. 



6 CHANGING THE BLACK HOLE MASS 

It is instructive to study the SED produced by the jet in FSRQs of 
different black hole masses, scaling the relevant quantities with the 
Schwarzschild radius and the Eddington luminosity. To this aim 
we show in Fig.QT]a sequence of SED with the black hole mass 
ranging from 3 x 10 7 to 3 x 10 9 Mq. We assume that -Rdiss is 
always at 500 Schwarzschild radii and that Ld = 0.067LEdd, with 
a corona with an X-ray luminosity equal to one third that of the 
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Table 1. List of parameters used to construct the SED shown in Figg. 6-13. Col. [1]: figure number where the model is shown; Col. [2]: dissipation radius in 
units of 10 15 cm and (in parenthesis) in units of Rs', Col. [3]: black hole mass in solar masses; Col. [4]: size of the BLR in units of 10 15 cm; Col. [5]: power 
injected in the blob calculated in the comoving frame, in units of 10 45 erg s — 1 ;Col. [6]: accretion disk luminosity in units of 10 45 ergs -1 and (in parenthesis) 
in units of I/Edd; Col. [7]: magnetic field in Gauss; Col. [8]: bulk Lorentz factor at -Rdiss; Col. [9]: viewing angle in degrees; Col. [10] and [11]: break and 
maximum random Lorentz factors of the injected electrons; Col. [12] and [13]: slopes of the injected electron distribution [Q(~f)] below and above 71,; Col. 
[14]: redshift; Col. [15]: some notes. For all cases the X-ray corona luminosity Lx = 0.3-Ld- Its spectral shape is assumed to be oc u~ 1 exp(— hu/150 keV). 



disk. We assume that the power in relativistic electrons injected 
into the dissipation region of the jet scales with the black hole mass 
as P( — 10 44 Af9 erg s _1 . Electrons are injected between 71 = 10 2 
and 72 = 10 4 in all cases. The bulk Lorentz factor is F — 15 for 
all cases and the redshift is 2 = 3. At the assumed i?di ss the jet has 
already reached its maximum F, and at these distances the Poynting 
flux is assumed to scale as Pb = 6.7 x 10 -3 Z/Edd- We then have 



Ub oc 



M 



1 



r 2 r 2 

1 -"-diss 



(36) 



which follows from the assumption of Rdi SB /Rs =const and 
r =const. For all our cases, the dissipation occurs within the BLR, 
which yields a constant radiation energy density (7 B lr since F is 
the same. As a consequence, the ratio Lec/ L ayn ~ <7blr/^b oc 
M. This is the reason of the increasing dominance of the inverse 
Compton emission increasing the black hole mass. This implies 
that blazars with large black hole masses should preferentially be 
more Compton dominated, and therefore more easily detected by 
the Fermi satellite. In fact, in Fig.QTJone can see the 5a sensitivity 



of Fermi for 1 year of operation (grey line), suggesting that, at high 
redshifts, the detected blazars will preferentially have large black 
hole masses. 



The importance of the EC relative to the SSC emission in- 
creases with the black hole mass, hardening the X-ray spectral 
shape. The SSC and EC components are shown separately for the 
SED corresponding to M = 3 x 1O 7 M0, to illustrate the im- 
portance of the SSC flux. For larger masses the SSC components 
becomes relatively less important than the EC one. For the SED 
with M = 3 x 10 9 A/q we show the effects of neglecting the 7- 
7 absorption and the consequent reprocessing (dashed line). One 
can see that the primary (i.e. neglecting pairs) spectrum as a rather 
sharp cut-off due to the Klein-Nishina limit given by Eq.|29l 
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Figure 11. The observed SED for different black hole masses (from 3 X 10 7 
to 3 X 10 9 A/q, as labelled) assuming that the dissipation takes place at 
500 i?s, that the accretion disk luminosity is = 0.067LEddi ar| d that 
the power injected in the jet dissipation region scales as P.' = 10 44 Afg 
erg s . The bulk Lorentz factor is kept fixed at V = 15. The particles 
are always injected between 71 = 10 2 and 72 = 10 4 . The grey line 
is the 5a detection sensitivity of Fermi, after 1 year of operation. For the 
M = 3 X 10 7 A/q case we show, besides the total spectrum, the SSC 
and the EC components separately. This illustrates the importance of the 
SSC process for low values of the black hole masses. The dashed line (for 
the case with M = 3 X 10® Mq) shows the spectrum neglecting photon- 
photon absorption and pair reprocessing. The received flux is calculated 
assuming that all sources are at z = 3. 



Figure 12. SED of 3C 279 in different states together with the correspond- 
ing models, whose input parameters are listed in Tab. [TJ See Ballo et al. 
(2002) and references therein for the sources of data points. Note that the 
high energy data-points (above 100 GeV) have been de-absorbed accord- 
ing to the Primack, Bullock & Somerville (2005) model used in Tavecchio 
& Mazin (2009). According to Eg. 1301 the flux at these energies can be pro- 
duced at relatively large distances from the black hole, beyond the BLR. 
In this case the main contribution to the seed photons is coming from the 
IR torus (solid light grey line, orange in the electronic version). In the low 
state (blue triangles) the SSC process is dominating the 2-10 keV X-ray 
flux. For this state only we show the SSC (long dashed line) and the EC 
flux (dot-dashed line) separately. 



7 SOME ILLUSTRATIVE EXAMPLES 
7.1 TeV FSRQs: the case of 3C 279 

3C 279 has been recently detected in the TeV band (Albert et al. 
2008), although its redshift, z = 0.536, implies a strong absorption 
of high energy 7-rays by the IR cosmic background. This demon- 
strates that also powerful blazars emit at large energies, up to the 
TeV band, even if the peak of their high energy hump may lye in the 
MeV-GeV band. Besides the consequences that this result has on 
the cosmic background, discussed in Sitarek & Bednarek (2008); 
Tavecchio & Mazin (2009); Liu, Bai & Ma (2008), we would like 
to discuss here another consequence, which concerns the primary 
spectrum of the source. 

We have discussed previously that if the bulk of the inverse 
Compton spectrum uses BLR photons as seeds, then we expect a 
steepening of the intrinsic spectrum at ~ 15/(1 + z) ~ 10 GeV 
due to Klein-Nishina effects. This is shown in Fig.[T2]as the darker 
line passing through an archival EGRET spectrum of the source 
(squares). For this model, in fact, the dissipation region is within 
the BLR, which is then giving most of the seed photons used for the 
scattering (see Tab. [TJ for the all the input parameters). The rather 
abrupt cut-off seen at ~ 2.5 x 10 24 Hz ~ 10 GeV is due this 
effect, and not to 7-7 internal absorption, which, in this particu- 
lar example, is unnoticeable. To produce photons of larger energies 
and evade the Klein-Nishina limit we are forced to locate i?di ss be- 
yond the BLR, and then use the IR photons produced by the torus. 



This is what the model does (lighter grey line, orange in the elec- 
tronic version). 

Note that the presence of a very high energy component in 3C 
279 was predicted (before detection) by Blazejowski et al. (2000), 
who included the presence of a IR emitting torus (see their Fig. 4). 
The main difference of our modelling with respect to Blazejowski 
et al. (2000) is in the assumed temperature of the torus, assumed 
to be larger in that paper, implying a smaller size and an enhanced 
U' m . 

If the jet of 3C 279 is "canonical" in the sense described in 
this paper, then a larger Ra ss means a smaller magnetic field, and 
a larger bulk Lorentz factor (if it has not yet reached its maximum 
value). Then, according to these ideas, we show the entire modelled 
SED to compare it with the "high EGRET state" (squares) and the 
"low EGRET state" (triangles). For the latter model we have as- 
sumed a smaller iidiss and P{. This implies a small bulk Lorentz 
factor (hence a decreased importance of the external seed photons) 
and a large magnetic field, resulting in a SED of equal synchrotron 
and inverse Compton power. 

3C 279 is one of the best studied 7-ray blazars, partly because 
it was very active during the observations of EGRET. It should 
not be taken, however, as the prototypical high power FSRQ, since 
its disk emission is very modest, as also directly suggested by its 
SED in the low state (Pian et al. 1999; see the triangles in Fig. 
1121 see Ballo et al. 2002 and references therein for the data): for 
our models we have assumed a black hole mass M = 1O 9 M 
and Ld = 0.02LEdd- It is therefore instructive to show example 
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Figure 13. Top panel: the SED of one of the most distant blazars, together 
to 3 different models with different i?<jiss ( see Tab. [T] for the set of pa- 
rameters) to illustrate possible different states of the source. See Fabian et 
al. (2001) and references therein for the sources of data, and Celotti et al. 
(2007) for further discussion about this source. 



of more powerful blazars, with accretion disk emitting close to the 
Eddington limit. 
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Figure 14. Top panel: the SED PKS 2149-307 together to 3 different mod- 
els (see Tab.[T]for the set of parameters) to illustrate possible different states 
of the source. These relatively distant, close to Eddington, high black hole 
mass blazars should be at the extreme of the blazar sequence, showing a 
high energy peak in the 100 keV-1 MeV energy band. Note that if the dis- 
sipation takes place very close to the black hole, when the jet is still accel- 
erating and with a strong magnetic field, the resulting spectrum becomes 
unconspicuous at high energies, even if the intrinsic dissipated power is the 
same of the higher states. 



7.2 MeV blazars at large redshifts 

GB 1428+4217, at z = 4.73, is the second most distant blazar 
known (the most distant is Q0906-6930, with z = 5.47). Its SED 
is shown in the top panel of Fig.[T3] together with the fitting model 
(black line). The figure shows that this source has a prominent UV- 
bump, with a luminosity around 10 47 erg s _1 (see also Tab.[TJ. 

According to our model, the peak of the high energy hump of 
the SED during the BeppoSAX observations (Fabian et al. 2001) 
is predicted to lye at ~ 1 MeV, locating this source at one ex- 
treme of the blazar sequence. The model assumes a black hole mass 
M = 1.5 X 10 9 M o , an accretion disk emitting at ~ 2/3 of the 
Eddington ratio (the other third is emitted by the X-ray corona). 
The dissipation region is within the BLR, at 1500 Schwarzschild 



radii from the centre. The fit requires B = 1.2 G, corresponding 
to P B = 5.4 x 10 45 erg s~\ and L d = 1.35 x 10 47 erg s _1 . 
The bottom panel of the same figures shows the radiation and mag- 
netic energy density profiles for the parameters used to model this 
blazars. The mid grey vertical line corresponds to the used value 
of -Rdiss- This figure shows that U^ LR is about two orders of mag- 
nitude larger than Ub, corresponding to the ratio of the external 
Compton to synchrotron luminosities. 

For illustration purposes, we also show the SED correspond- 
ing to increase (decrease) iidiss tenfold (see the vertical grey lines 
in the bottom panel, and the corresponding SED in the top panel). 
Shifting R diBB to 1.5 x 10 4 i?s the high energy peak moves to lower 
frequencies, as the relevant seed photons are softer, and at the same 
time the X-ray spectrum hardens. This is obtained even if the par- 
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tide injection function is unchanged, and is due to the incomplete 
cooling occurring for the larger i?diss case. Since 7 coo i shifts from 
1 to ~34, in the X-ray band we see the emission from the un- 
cooled electron populations, that retains the original, "injection", 
slope s = 0. Note that the high energy peak is now slightly be- 
low 100 keV, and this would make the source even more extreme 
in terms of the blazar sequence. The EC to synchrotron luminosity 
ratio is still ~ 100, corresponding now to the ratio between f/ IR 
and t/g, as can be seen in the bottom panel. 

Instead, if i?di ss is at 150 Schwarzschild radii, the SED 
changes more dramatically. Since in this case the bulk Lorentz fac- 
tor is smaller and the magnetic field larger, the inverse Compton 
and synchrotron powers become comparable (see the bottom panel, 
showing that for this i?di ss we have U-g ~ U' d ). The X-ray spec- 
trum brightens and softens considerably, in the IR-UV band the 
flux increases (because of the increased magnetic field) even if the 
bolometric observed luminosity decreases because of the decreased 
Doppler boosting. The fact that we did not see this kind of SED in 
GB 1428+4217 suggests that this state rarely occurs. This is partic- 
ularly true considering that this blazars was not discovered because 
it was particularly bright in hard X-rays or in the 7-ray band, so 
there was no bias against a soft X-ray spectrum. 

We lastly consider PKS 2149-307 at z = 2.345, as observed 
by the XRT and BAT instruments onboard Swift during its first 
9 months of observations (Sambruna et al. 2007). Fig. 1141 shows 
its SED (top panel) and the profiles of the radiation and mag- 
netic energy densities for the considered models (bottom panels). 
As done for GB 1428+4217, we show what we considered the 
best fitting model (black line in the top panel) corresponding to 
.Rdiss = 9.6 x 10 17 cm, (corresponding to 800 Rs for the assumed 
black hole mass of M = 4 x 1O 9 M0), and also the models corre- 
sponding to a tenfold increase (decrease) of i?di ss (light grey lines 
in the top panel and vertical lines in the bottom panel). The BAT 
data have large error bars, precluding the possibility to firmly claim 
that its high energy peak is within the BAT energy range (namely, at 
~ 100 keV), but this possibility is indeed suggested by the present 
data. The SED of this source can be explained by a set of param- 
eters similar to the ones chosen for GB 1428+4217, and similar 
considerations apply. Increasing i?di ss tenfold makes the IR radia- 
tion from the torus to dominate, leaving almost unchanged the ra- 
tio between the radiation and the magnetic energy densities (and 
thus the corresponding EC to synchrotron luminosity ratio). Again 
a decrease in i?di ss by a factor 10 has a dramatic impact on the pre- 
dicted SED, dominated in this case by the synchrotron luminosity 
(see also the bottom panel). However, having only data up to the 
medium-soft X-ray energy range, it would be almost impossible 
to tell what is the real appearance of the bolometric SED of the 
source, as almost all the changes occur above 10 keV. 



8 CANONICAL JETS AND THE BLAZAR SEQUENCE 

We can ask if our "canonical" jet can reproduce the phenomeno- 
logical blazar sequence as proposed by Fossati et al. (1998) (see 
also Donato et al. 2001 and Ghisellini et al. 1998). This sequence 
was constructed by taking flux limited samples of blazars (in the 
radio and X-ray bands), dividing the sources in radio luminosity 
bins, and averaging the flux of the sources in each (radio lumi- 
nosity) bin at selected frequencies. The resulting data are shown 
in Fig. [15] As discussed in Ghisellini & Tavecchio (2008) and in 
Maraschi et al. (2008), this sequence does not pretend to describe 
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Figure 15. The blazar sequence (Fossati et al. 1998; Donato et al. 2001) 
interpreted in the framework of our canonical jet scenario. Parameters are 
listed in Tab.[T] The SED changes according to changing the accretion rate 
and the power of the jet, and assuming that below some critical accretion 
rate the accretion regime changes regime, becoming very radiatively inef- 
ficient. In out case, this occurs for the two least powerful SED, that should 
corresponds to low power, line-less BL Lacs. 



"the average blazar", since it likely represents the most beamed 
sources, successfully entering in the flux limited sample they be- 
long to. Nevertheless it is appropriate to ask if our "canonical" jet 
can reproduce this sequence without strong modifications of the 
basic assumptions. Fig.[T5]then shows the SED resulting from our 
modelling, whose parameters are listed in Tab. Q] As can be seen 
the agreement is quite good, and it is achieved by assuming (as in 
Ghisellini & Tavecchio 2008), that the accretion disk becomes ra- 
diatively inefficient below some critical accretion rate (thus around 
a critical luminosity in units of the Eddington one, that we take as 
a few x 10 -3 ). The first three, more powerful, SED correspond to 
sources with a standard disk, while for the two SED at lower lumi- 
nosities we have "switched-off" all external radiation. The black 
hole mass is 3 x 1O 9 M0 for the two most powerful SED, and is 
10 9 Mq for the other three. The luminosity, in the comoving frame, 
injected in relativistic electrons is monotonically decreasing, as the 
value of the magnetic field. The dissipation radius i?di ss , in units 
of the Schwarzschild radius, changes slightly, but less than a fac- 
tor 2.5. The energy of the injected particles need to increase as the 
cooling decreases. The decreased cooling also makes 7 coo i to in- 
crease. These effects are very important for the two "BL Lac" SED, 
and much less for the other three most luminous SEDs. 

We can conclude that the blazar sequence, as illustrated in Fig. 
1151 can be interpreted as a sequence of jet powers, in sources having 
more or less the same (large) black hole mass. The jet power cor- 
relates with the accretion disk luminosity and below some critical 
value the disk changes accretion regime, becoming radiatively inef- 
ficient and unable to photo-ionise the broad line clouds. A intrigu- 
ing manifestation of this behavior is shown by the bright blazars 
recently detected by the Fermi satellite (Abdo et al. 2009) dis- 
cussed in Ghisellini, Maraschi & Tavecchio (2009), and mirrors 
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what has been suggested to occur for the FRI-FRII radio-galaxy 
divide (Ghisellini & Celotti 2001). 



9 CAVEATS 

Our aim was to describe the main properties of the produced spec- 
trum in a "canonical" jet. Specific sources can of course deviate 
somewhat by our description. For instance: 

1/2 

• Assuming -Rblr oc L d , may be approximately obeyed 
on average, but specific sources might behave differently. The as- 
sumed jet geometric profile and the assumption that the jet, at its 
start, is magnetically dominated could be two oversimplifications. 
We hope to test these assumptions by direct observations, since the 
combined efforts of the Fermi and Swift satellites can give us really 
simultaneous spectra on a large frequency band. 

• We have assumed that the "dissipation" mechanism, convert- 
ing some fraction of the jet power into radiation, has the form of an 
injection of primary leptons lasting for a light crossing time. With 
this assumption, we have derived the particle distribution at this 
time, which corresponds to the maximum produced flux. However, 
the injection can last for a longer or shorter time. Our framework is 
appropriate for describing a "snapshot" of the SED, and not for a 
time-dependent analysis. 

• We assume that the magnetic field in the dissipation region is 
the same as the magnetic field transported by the jet. It is not am- 
plified by e.g. shocks, nor it is reduced by e.g. reconnection events 
(see e.g. Giannios, Uzdensky & Begelman 2009). Although we can 
assume any value of the magnetic field in the dissipation region (i.e. 
we can fix any value in the numerical code calculating the SED), 
we prefer this choice for all the examples shown, because of sim- 
plicity. 

• We focused on high power blazars, thought to have "standard" 
accretion disks and broad line regions, and therefore we did not 
discuss in detail the expected SED from low power (and TeV emit- 
ting) BL Lacs. However, we could simply extend our study to these 
sources if we "switch-off the external radiation components, as 
expected if the accretion disk becomes radiatively inefficient, mak- 
ing very weak or absent broad emission lines. In Fig.[T5]we indeed 
show two examples of the expected spectra in this case. 

• All the shown SED do not take into account the absorption 
of 7-rays due to the IR, optical and UV background, but only the 
"internal" absorption caused by the radiation fields existing close 
to the emitting blob. 

• We have always assumed a relatively small viewing angle 
(9 V ~ 3°), of the same order of 1/r. Since Y and 9 V are two sepa- 
rated input parameters, we could change 9 V (and thus the Doppler 
factor 5) for a fixed T. However, this would introduce a delicate is- 
sue, since the external radiation field, in the comoving frame, is not 
isotropic. If the emitting blob is inside the BLR, for instance, then 
all external photons appear to come from one direction, and we 
could apply the formalism of Dermer (1995) to describe the pat- 
tern of the scattered radiation. Note that, in this case, if 6 V — 1 /Y 
(and thus 8 = T) one has the observed luminhosity as for radiation 
isotropic in the rest frame. On the other hand, if the emitting region 
is beyond the BLR, then the arrival directions of the BLR exter- 
nal photons are spread, and it is not trivial to reconstruct the exact 
pattern. But in this cases the BLR component is hardly dominant 
with respect to the IR radiation from the torus. We have decided, 
for simplicity, to use 9 V ~ 1 /Y, and to treat the emitted radiation 
pattern is the same way as for an isotropic seed photon distribution. 



• We neglected, for simplicity, the possibility of a jet composed 
by a fast spine surrounded by a slower layer. This structured jet 
have been proposed for low power BL Lacs and FR I radio-galaxies 
(Ghisellini, Tavecchio & Chiaberge 2005), and can ease some prob- 
lems in explaining the SED properties in these sources. While in 
high power blazars there is no compelling need, yet, for such a 
structure, its inclusion in the present study would make our descrip- 
tion much more complex and model-dependent, given the freedom 
to choice the parameters for the layer emission. 

• For simplicity, we neglected the possibility that the BLR and 
the associated inter-cloud material can absorb and re-emit, or scat- 
ter, part of the synchrotron radiation produced by the jet (i.e. the 
"mirror model", see Ghisellini & Madau 1996, and some criticism 
by Bednarek 1998; Bottcher & Dermer 1998). 

• The jet, even if cold (i.e. before dissipation), is relativistic, 
and the carried leptons could scatter the external radiation by bulk 
Comptonization. This process, proposed by Begelman & Sikora 
(1987) and Sikora et al. (1994) and studied in detail in Celotti, Ghis- 
ellini & Fabian (2007), could produce a black-body like component 
in the X-ray band, whose level depends on how the jet accelerates 
and on the amount of leptons carried by the jet. We have neglected 
this component for simplicity, although it can be a very important 
diagnostic for deriving, at the same time, the bulk Lorentz factor of 
the jet and its matter content. 



10 CONCLUSIONS 

The relatively recent observations of blazars, especially in the y— 
ray band, have disclosed some of the crucial ingredient for our un- 
derstanding of jets, namely the bolometric power output and the 
real shape of the produced SED. Although there is still discussion 
about the origin of the high energy emission of blazars, the rel- 
atively simple one-zone leptonic model was rather successful to 
interpret the existing data of specific sources. 

In this paper, rather than modelling single sources, we tried 
to propose a more general description of the emission produced by 
high power jets. For this aim we have first summarised in a simple 
way the expected radiation fields as a function of the distance from 
the black hole, and then we have studied how the expected SED 
changes if the main dissipation occurs at different locations in the 
jet. We have done this by assuming rather reasonable prescription 
for the accretion disk and its X-ray corona emission, and using 
reasonable prescriptions for the radiation fields corresponding to 
the broad line region and the infrared torus. We have also assumed 
that the jet conserves its Poynting flux at all distances, thus fixing 
the profile of the magnetic field. Since all these prescriptions cor- 
respond to rather uncontroversial assumptions, we believe that they 
should well characterise the properties of an average, canonical, 
large power jet. 

The main conclusions of our work are: 

• The magnetic energy density is always smaller than the exter- 
nal radiation energy densities, except at the very beginning of the 
jet and at the kpc scale. Beyond these scales, the contribution of the 
cosmic microwave background becomes dominant. 

• The inclusion of the X-ray radiation from the disk corona al- 
lowed to quantify the effect of pair production and reprocessing 
for small i?di ss - This is found to be severe if the primary spec- 
trum emits most of its power at energy above the pair production 
threshold. This requires that the bulk Lorentz factor is large even 
in the very vicinity of the black hole. In this case there is a trans- 
fer of power from the 7-ray to the X-ray band. Our quantitative 
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analysis confirms earlier, more qualitative, statements concerning 
the non-dissipative nature of the inner (i.e. up to a few hundreds 
Schwarzschild radii) jet. 

• If instead the jet accelerates gradually (i.e. as T k R}l 2 up 
to a maximum value), the magnetic energy density dominates the 
cooling at the start of the jet, resulting in a synchrotron dominated 
SED, with a very small fraction of power being emitted above the 
pair production energy threshold. These SED would hardly be de- 
tectable, since the small Y means a much reduced Doppler enhance- 
ment of the flux, and thus these components are easily overtaken by 
the emission at larger distances (that have a larger T). In any case, 
early dissipation (i.e. small -Rdiss) is not an efficient process to cre- 
ate electron-positron pairs. 

• Besides the inner part of the jet, the "internal" pair absorption 
(i.e. calculated within the emitting blob) can be particularly impor- 
tant also if dissipation takes place close to (but within) the BLR and 
close to (but within) i?diss ~ Rm- In these cases, however, a small 
fraction of the emitted power can be absorbed, and the reprocessing 
is modest. 

• The Klein-Nishina decrease (with energy) of the scattering 
cross section imprints a characteristic steepening in the spectrum 
above ~ 10 GeV if the dissipation takes place within the BLR, and 
at ~ 1 TeV if the dissipation occurs beyond the BLR but within 
Rm- 

• The above point has important consequences on the possible 
use of the high energy (GeV) data of blazars to study the optical- 
UV background. In fact the signature of the 7-7 absorption due to 
this component would be a rather sharp steepening of the received 
spectrum, that would be very similar to the intrinsic steepening of 
the primary spectrum due to Klein-Nishina effects ;/ most of the 
dissipation takes place within the BLR. In this case also the "inter- 
nal" 7-7 absorption due to BLR photons will be relevant. It will be 
very difficult to disentangle these effects. On the other hand, if the 
observed spectrum continues unbroken above 15(1 + z) GeV, we 
can infer that most of the dissipation has occurred beyond -Rblr- 
In this case both Klein-Nishina and "internal" absorption effects 
are unimportant up to the TeV band. These therefore would be the 
best candidates to study the opt-UV cosmic backgrounds. 

• The magnetic field decreases with distance, while the C/ B lr 
is approximately constant up to -Rblr,, and U( R is approximately 
constant up to Rm. At these distances (i.e. Rdi ss ~ -Rblr and 
-Rdiss ~ Rm) the Compton dominance is very large. Existing ob- 
servations of large power blazars already indicate that these blazars 
are characterised by a weak synchrotron component, a strong ther- 
mal, disk component (unhidden by the synchrotron flux) and a very 
large high energy flux. These are the sources with the largest Comp- 
ton dominance, strongly suggesting that -Rdiss is relatively large 
(i.e. close to the -Rblr or Rm). This helps to understand why in 
these sources the soft X-rays are relatively weak and very hard, de- 
manding a very weak SSC flux: since the magnetic field is small, 
not only the synchrotron, but also the SSC flux is small with respect 
to the EC radiation. 

• Concerning the above point, the fact that the high power 
blazars we already know of have the mentioned properties does 
not imply that they preferentially dissipate at large distances with 
respect to less powerful sources, since we may have selected them 
on the basis of their hard X-ray or 7-ray flux (therefore taking ad- 
vantage of the large Compton dominance they have when -Rdiss is 
large). In other words: they may dissipate at different -Rdiss at dif- 
ferent times, but we select them only when -Rdiss is large. 

• Although most of the dissipation should take place in one zone 
of the jet (with .Rdiss -500-1000 Rs), it is very likely that dissi- 



pation occurs also at larger distances, albeit with less power. The 
flux originating in these zones can be important in some frequency 
bands, as the soft X-ray one, with important consequences on the 
observed variability. Instead, in the 7-ray energy band (at least 
when the 7-ray dominates the bolometric output) most of the flux 
should originate in one zone. 

• Blazars with large masses should be characterised by the 
largest inverse Compton to synchrotron ratios, and thus they will 
be well represented in 7-ray all sky surveys, as performed by the 
Fermi satellite. 

• The most powerful blazars, with the most powerful associated 
accretion disk, should have a high energy hump peaking in the 100 
keV-1 MeV energy range. The exact location of the peak frequency 
depends somewhat on -Rdiss, being smaller if dissipation occurs be- 
yond the BLR. These blazars should form the high power extension 
of the blazar sequence. They are not easily detectable by the Fermi 
satellite, since their MeV peak implies relatively small fluxes in the 
GeV band. The brightest of them have already been revealed by the 
serendipitous survey of the BAT instrument onboard Swift. 

• The phenomenological blazar sequence as presented in Fossati 
et al. (1998) can be interpreted as a sequence of "canonical jets" in 
sources with more or less the same (large) black hole mass, but with 
different accretion rates and different jet powers. 
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